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The observed steep kinetic scale turbulence spectrum in the solar wind raises the question of how that 
turbulence originates. Observations of keV energetic electrons during solar quiet time suggest them as a 
possible source of free energy to drive kinetic turbulence. Using particle-in-cell simulations, we explore 
how the free energy released by an electron two- stream instability drives Weibel-like electromagnetic 
waves that excite wave-wave interactions. Consequently, both kinetic Alfvenic and whistler turbulence are 
excited that evolve through inverse and forward magnetic energy cascades. 

DOI: 10.1 103/PhysRevLett.l 12.061 101 PACS numbers: 96.60.Vg, 52.25.Dg, 52.35.Ra, 94.05.Lk 


The observations of solar wind turbulence have shown 
that at scales approaching the ion kinetic scale where the 
ions and electrons are decoupled and the kinetic effects 
must be considered, the power spectrum of magnetic 
fluctuations, which in the inertial range follows the 
Kolmogorov scaling cx kr 5 / 3 , is replaced by a steeper 
[1-3] anisotropic scaling law B 2 ^ oc k~]_ a , where a is a 
number larger than 5/3. It is found that the observed 
spectral index is a ~ 2.7 but not universal and varies in a 
range [2,4]. Magnetic fluctuations with about a tenth of ion 
gyrofrequency propagating nearly perpendicularly to the 
solar wind magnetic field are identified as kinetic Alfvenic 
waves (KAWs) [4-9], and the break frequencies of the 
magnetic power spectra appear to follow the ion inertial 
length [3,10,11]. The origin of the KAW turbulence is still 
unknown. In this Letter, we address the origin of kinetic 
turbulence by proposing a source of free energy that has not 
been explored previously. For the first time, we find that an 
inverse energy cascade appears to play a crucial role in 
generating both KAW and whistler turbulence. 

Observations using the STEREO spacecraft have found 
that, even during quiet-time periods, the solar wind con- 
tains a previously unknown electron population different 
from the core solar wind, called “superhalo electron,” with 
energy ranging from to 20 keV [12,13]. One possible 
origin of the superhalo electrons is the escaping nonthermal 
electrons related to coronal nanoflares in the quiet solar 
atmosphere [12,14]. The relative drift of these nonthermal 
electrons to the background solar wind can drive an 
electron two-stream instability in a neutral current [15] 
and release the free energy to the solar wind. The impact of 
this unstable process on the solar wind has so far not been 
studied. In this Letter, using particle-in-cell simulations, we 
investigate how the rapidly released energy drives Weibel- 
like electromagnetic waves. The wave-wave interactions on 
ion inertial scales d t = c/co pi and electron inertial scales 
d e = c / co pe generate KAW and whistler turbulence 
through both forward and inverse energy cascades. At 


the end of this Letter, we will compare the testable features 
produced by this model with observations. 

We initialize the 2.5 D particle-in-cell simulations in the 
solar wind frame of reference with a uniform magnetic field 
B = B 0 x. Both the ion and the electron densities are 
uniform. The initial ion velocity distribution function is 
a single isotropic Maxwellian. The electron velocity dis- 
tribution function is a core-beam isotropic bi-Maxwellian. 
The core is the solar wind electrons and the beam is the 
energetic electrons. Their relative drift is along B 0 : 
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where v\ = v 2 ey + v 2 ez , 8 = n b /n 0 . n 0 is the solar wind 
density and the density normalization unit, n b is the density 
of beam electrons. T c is the temperature of the core, T b is 
the temperature of the beam, v cd is the drift of the core, and 
v bd is the drift of the beam. The drift velocities satisfy 
(1 — 8)v cd = —8v bd to maintain null current. v bd = 
I2v te = 60 v A , where v A is the Alfven speed and v te = 
VkTJ m e is the thermal velocity of the core electrons. The 
energy of these beam electrons will be released and join the 
core electrons at energy ~kT c . We choose 8 = 0.1, because 
at ~10 keV, or ~10 3 kT c , the superhalo electrons have a 
density of ~10 -6 of the solar wind density [13], and we 
assume the kinetic energy flux density of beam n b v\J2 is 
constant. The speed of light in these simulations is chosen 
to be c — 100i; A and the mass ratio is m i lm e = 100. 
The ion temperature T l — T c . The boundaries are periodic 
in both directions with a box size L x = L y = 32 d t . 
The total number of cells in each dimension is 10 240 
and the total number of particles is ~10 10 . The total 
simulation time co pe t = 10560. The electric field is nor- 
malized by E 0 = v a B 0 /c. We take kT b — 2kT c = 0.5m^ A 
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FIG. 1 (color online), (a) Time evolution of energy of (8 El) 
(black line), (<55 2 ) (red line), (8 By) (green line), and (8B\) (blue 
line). The embedded plot is an expanded view of the time 
evolution from co pe t — 0-230. The orange line indicates 
co pe t = 230. (b) Power spectrum of \8E x (k x , k y )\ 2 at co pe t = 
24 on a logarithmic scale. 


and j5 = kT c /B\ = 0.25 estimated from the solar wind ft 
observations at 0.3 AU [11]. 

Electron two-stream instability occurs early at 
c o pe t = 24, as shown in Fig. 1(a), and (8 El) (solid black 
line) quickly reaches a peak at co pe t « 50, where () denotes 
the average over xy. At co pe t = 200, the drift of the beams 
decreases from 60v A to ~20v A , and 8E X decreases by 
nearly a factor of 20 and then stays nearly constant. 
The growth rate of the electron two-stream instability at 
co pe t = 24 is close to the cold plasma limit of 
Yb ~ { n b/^ n oY^ 3(D pe ~ 0Aco pe . The fastest growing mode 
kf x = co pe /v db ~ 11/ dj is consistent with the spectrum of 
1 8E x (k x ,k y )\ 2 at co pe t = 24, as shown in Fig. 1(b). 

The fast growth of 8E X generates an inductive magnetic 
field B z that satisfies B z ~ ( 8E x Ay/cAt ) ~ 0.245 0 , which 
is close to the middle value of B z shown in Fig. 2(a), where 
we take Ay ~ Af tX = 2n/kf x ~ 3d e , At ~ l/y b ~ 2.5co p l. 
The middle value of 8E X ~ 20 during the instability is 
estimated from Fig. 1. The internal energy density re- 
leased per wavelength is ~m e n b (Av db ) 2 Ay / (2At)~ 
O.2n 0 m i v\co pe 2j\, where Av db ~60v A . Around 10% is con- 
verted into magnetic energy B 1 /%n ^ ^S)^n^miv\(D pe X~j\ 
at the end of the two-stream instability, while nearly 90% is 
converted into the thermal motion of trapped electrons [16]. 

The electric current density j ex produced by the induc- 
tive magnetic field becomes as important as the displace- 
ment current when the two-stream instability starts to 
decay. Then j ex drives a Weibel-like instability that gen- 
erates nearly nonpropagating transverse electromagnetic 
waves. The variances ((55 2 ) 1 / 2 and (<55 2 ) ly/2 in Fig. 1(a) 
reach a second peak at co pe t^612. The second peak is 
much higher than the first peak produced by the two-stream 
instability. The variance (<55 2 ) 1//2 follows ((55 2 ) 1 / 2 closely, 
while the variance ((55 2 ) 1 / 2 reaches its peak at a slightly 
later time. A significant change from electrostatic waves to 
transverse electromagnetic waves can be seen in the 
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FIG. 2 (color online). Images of B z /B 0 at co pe t — 24 (a), 480 
(b), 2424 (c), and 10 560 (d). Refer to the movie in the 
Supplemental Material [24]. 



evolution of j ex , shown in Fig. 3. The j ex wave vector 
induced by the two-stream instability is along v. Gradually, 
the wave vector rotates so that it is parallel to y, which 
indicates the generation of electromagnetic fluctuations in 
B z that align along y. The wavelength of B z fluctuations 
increases to half d h as seen in Fig. 2(b) at co pe t = 480, near 
the peak of the Weibel-like instability. 

The decay of the Weibel-like instability enhances the 
interactions between the localized currents and the nearly 
nonpropagating transverse electromagnetic waves. This 
process breaks up the transverse waves and produces 



x / d j 


FIG. 3 (color online). The transition of j ex wave patten when 
Weibel-like instability occurs, (a) j ex at co pe t — 96, the late stage 
of two-stream instability, (b) j ex at co pe t — 168, the transition 
stage from the two-stream instability to the Weibel-like insta- 
bility. (c) j ex at co pe t = 240, the beginning of the Weibel-like 
instability. 


061101-2 







PRL 112, 061101 (2014) 


PHYSICAL REVIEW LETTERS 


week ending 
14 FEBRUARY 2014 


randomly propagating waves as shown in Fig 2(c). From 
c o pe t = 2400, the wave-wave interactions dominate the 
dynamics. The wave-wave interactions lead to a momen- 
tum transfer from perpendicular to parallel magnetic field. 
As a result, parallel propagating waves appear, which is 
consistent with the fact that a peak appears in ( SB 2 ) 1 ^ 2 at 
co pe t = 2400 (Fig. 1). Finally, at co pe t = 10560, nearly 
perpendicular propagating waves with angle > 89° and 
nearly parallel waves are excited [Fig. 2(d)]. 

The wave-wave interactions drive a bidirectional energy 
cascade. The perpendicular magnetic wave energy is now 
transferred from the electron inertial scale back to the ion 
inertial scale, and the parallel magnetic wave energy is 
transferred from the ion inertial scale down to the electron 
inertial scale. The 2D power spectra of 5B Z at co pe t = 24, 
480, 2424, and 10 560 are shown in Figs. 4(a)-(d), 
respectively. At co pe t = 24, we see only a transverse mode 
peaked k y d t ~ 10, i.e., k y d e ~ 1, which is consistent with 
the wavelength of the inductive magnetic field B z that was 
produced by the two-stream instability. At oo pe t = 48, 
the Weibel-like instability generates a transverse electro- 
magnetic magnetic field with longer wavelengths. At 
oo pe t > 2424, wave-wave interactions occur in which a 
parallel branch is produced while the wave number of the 
perpendicular branch decreases. At the end of the simu- 
lation, magnetic powers are concentrated in two branches in 
the energy spectrum: the nearly perpendicular branch with 
k x d t < 1, and the parallel branch with k y d t < 2. We study 
the time evolution of the magnetic components of waves; 
the results show that both wave types are right-hand 
polarized. During the evolution, the magnetic wave-wave 



FIG. 4 (color online). Power spectra \B z (k x , k y )\ 2 on logarith- 
mic scale at co pe t — 24 (a), 480 (b), 2424 (c), and 10 560 (d). 


interactions form localized thin current sheets with widths 
from several d e to d h some of which might be caused by 
magnetic reconnections (see Fig. 2 in the Supplemental 
Material [24]). 

The frequency of the nearly perpendicular wave is 
around 0.2-0.3£^, where is the ion cyclotron frequency. 
From the dispersion relation of KAW given by the two- 
fluid equation [17] 


co 2 _ 1 + k 2 p 2 
k 2 x v\ 1 + k 2 <fi e ’ 

where p 2 s = d 2 v 2 J v\, we estimate k y d t < 8 for k x d t ~ 0.01 
and the electron thermal velocity is larger than the initial 
velocity v 2 e > T c /m e = 25 v\. The resulting KAW k y d t is 
consistent with the spectrum shown in Fig. 4(d). The 
frequency of the parallel branch is co ~ 10 and the wave 
number is kd t ~ 3 at co pe t = 10560, which satisfies the 
whistler wave dispersion relation co/£li = v A (kdi) 2 cos 6 
[18] for #~0. But at the transition time co pe t = 2424, 
kd t ~ 4, then 6 ~ 45° and k x d t ~ 2.5. Thus, the oblique 
whistler wave evolves to parallel. The ratio of 5B 2 /5B 2 ~ 1 
at the late stage (Fig. 1) implies that the turbulent magnetic 
energy is nearly equally distributed between KAW and 
whistler wave turbulence. 

Three- wave interaction ± k 2 = k 3 is the dominant 
process in wave-wave interactions and leads to the simul- 
taneous generation of KAWs and whistler waves. For 
perpendicular interactions, the major contribution is from 
Jfc^ w ± = kf™, where |^ h f ler | <<c |Jfc^ w |; thus, 

^jpw moves smaller wave numbers and the magnetic 
energy transfers from small scale to large scale. For parallel 
interactions, the major contribution is from fcj^ w zb 

^whisller = yfcwhistler ? where |^KAW| ^ ^ [whistler. thus? 

^whistler moves to larger wave numbers and the magnetic 

energy cascades down to the small scales. 

In Fig. 5, we show 1 D power spectra of the magnetic 
energy 5B 2 (k) vs k x (parallel spectrum) and 8B 2 (k) vs k y 
(perpendicular spectrum) at co pe t = 10560. Whistler wave 
energy cascades from ion to electron scales, and it is clear 
that the contribution to the parallel spectrum in k x on the 
ion scale is from whistler waves. The perpendicular 
spectrum has a bump at k y d t ~ 4-7 corresponding to 
the relic of magnetic energy injection by the two-stream 
instability at co pe t = 24 [Fig. 4(a)]. Then the wave-wave 
interactions inversely transfer the KAW energy to ion 
scale smaller than d t and generate the whistler waves 
[Fig. 4(d)]. Thus, both the whistler waves and KAWs 
contribute to the perpendicular spectrum on ion scale 
while only KAWs contribute to the perpendicular spec- 
trum on electron scale L e with d e > L e > p e . The spec- 
trum is much steeper on L e since the wave-particle 
interactions are much stronger. The parallel spectrum 
on a scale smaller than d e and the perpendicular spectrum 
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on a scale smaller than p e suggest exponential decays that 
imply that the dissipation processes are less space and 
time correlated. The plateaus between the power law 
and the exponential decays indicate that the energy is 
accumulated by the strong thermalization. 

After co pe t = 10560, the free energy is almost fully 
released and the induced turbulent scattering produces a 
nearly isotropic electron halo superposed over the core 
electrons. The energy exchange between particles and 
waves reaches a balance. The turbulence reaches its new 
steady state with P 2 + B 2 /Sn = const, where P is the total 
pressure of ions and electrons. The ratio of the amplitude of 
the magnetic fluctuations and the background magnetic 
field is about 0.2 and matches the current observations of 
solar wind kinetic turbulence. The decay rate of fluctua- 
tions is <$C 10 -4 £2 j in the last 3000 co~l ~ 2 O^ 1 , estimated 
from the simulation. This suggests that the kinetic turbu- 
lence will be preserved for a long time. If superhalo 
electrons are produced in the Sun, then kinetic turbulence 
is produced within a few solar radii. The resulting turbu- 
lence should be stable enough to travel to 1 AU based on 
the decay rate estimated from our simulation if we 
take £li ~ 1 Hz. 

Our model can naturally explain some important solar 
wind turbulence observations. (1) Current observations can 
be compared to the perpendicular power spectrum in Fig. 5 
on scales above d e . The spectral index of —2.2 in our 
spectrum agrees with those found in observations. In our 
simulation the KAW as well as the parallel propagating 
high frequency whistler waves contribute to the power 
spectrum. In Ref. [19], the authors suggest that ~10% of 
the solar wind data they analyzed consist of parallel 
propagating whistler waves as determined by their right- 
handed polarization, but more advanced observations are 
needed. (2) Since the growth rate of two-stream instability 
is related to co pe , the spectral breakpoints of KAWs follow 



FIG. 5 (color online). The ID spectra of 8B 2 {k) vs k x d t and 
k y d t . The blue short-dashed line is the k y space range for 
magnetic energy injection. In the simulation, we have d t — 2 p t 
and d e = 2 p e , and p t and p e correspond to k xy d t = 2 and 20 
respectively. 


the inertial length. This agrees with observations of KAW 
turbulence [3,10,1 1]. (3) At the final stage, enhanced by the 
relic parallel electric field from the two- stream instability, 
( | £|| | / \E ± | ) ~ 2-3, consistent with the observations that the 
parallel electric field is larger than the perpendicular 
electric field expected for KAW [20]. (4) During the 
evolution of the turbulence, microscopic current sheets 
with widths varying from several d e to d t are produced 
(Fig. 2 in the Supplemental Material [24]), consistent with 
the observations of kinetic scale current sheets discovered 
in solar wind turbulence [21,22]. (5) Our simulations show 
that a nearly isotropic halo is produced at the finale stage. 
Such halos are observed from 0.3 to 1 AU in slow wind [23] 
(Fig. 1 in the Supplemental Material [24]). The formation 
of a superhalo requires that the electron beam energy 
extend by more than an order of magnitude higher, 
rendering computations rather expensive due to the higher 
c/v A ratio and higher temporal and spatial resolutions 
required. 

The observations of kinetic turbulence on the electron 
scale would be more challenging since the power spectra on 
the electron scale is much steeper and quickly become 
exponential. The ongoing Magnetospheric Multiscale 
Mission might be able to detect the kinetic process on 
the electron scale at 1 AU. 

It is important to know how other possible turbulent 
processes in the solar wind affect the KAWs and whistler 
waves when they travel to 1 AU. The index of the power 
spectra might be affected more while other features 
produced in our model may be slightly or unaffected: 
the frequency breakpoints determined by ion inertial 
length, the enhanced parallel electric field, and the electron 
halo. The current observations of solar wind can reach 
0.3 AU. The near-future space missions Solar Probe Plus 
and Solar Orbiter can reach 10 solar radii and hence provide 
more rigorous constraints on this model. 

This model proposed is motivated by the observations of 
a superhalo. It is encouraging that the model could 
potentially link the existing observations of solar wind 
kinetic turbulence, the halo formation, and the electron 
acceleration and heating processes in solar corona into a 
coherent picture. More advanced studies will be carried out 
in the near future. 
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